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Abstract 

Bekenstein's model of a scalar field, 0, that affects the electromagnetic permeability 
(usually identified with "changing a" ) predicts tiny variations of the effective fine structure 
constant up to very high red-shifts, \a(z = 3.5)/ao — 1| < 10~ 10 , when the constraints from 
Eotvos-Dicke-Braginsky types of experiments are imposed. We generalize this model by 
allowing additional couplings of 4> to both a dark matter candidate and to the cosmological 
constant. We show that in a supersymmetric generalization of Bekenstein's model, the 
coupling to the LSP, which is assumed to contribute significantly to the dark matter density, 
can be up to six orders of magnitude stronger than the coupling to the baryon energy density. 
This allows one to evade the present limits on the non-universality of the gravitational 
attraction due to ^-exchange and at the same time accommodate the effective shift in a 
at the level of a(z = 3.5)/ a,® — 1 ~ 10~ 5 , reported recently from observations of quasar 
absorption spectra. 



1 Introduction 



Speculations that fundamental constants may vary in time and/or space go back to the 
original idea of Dirac [3|. Despite the reputable origin, this idea has not received much 
attention during the last fifty years for the two following reasons. First, there exist various 
sensitive experimental checks that coupling constants do not change (See, e.g. @). Second, 
for a long time there has not been any credible theoretical framework which would predict 
such changes. 

Our theoretical mindset, however, has changed since the advent of the string theory. 
One of the most interesting low-energy features of string theory is the possible presence 
of a massless scalar particle, the dilaton, whose vacuum expectation value defines the size 
of the effective gauge coupling constants. A change in the dilaton v.e.v. induces a change 
in the fine structure constant as well as the other gauge and Yukawa couplings. The 
stabilization of the dilaton v.e.v., which usually renders the dilaton massive, represents one 
of the fundamental challenges to be addressed before string theory can aspire to describe 
the observable world. Besides the dilaton, string theory often predicts the presence of other 
massless or nearly massless moduli fields, whose existence may influence particle physics 
and cosmology and may also change the effective values of the coupling constants as well. 

Independent of the framework of string theory, Bekenstein H formulated a dynamical 
model of "changing a". The model consists of a massless scalar field which has a linear 
coupling to the F 2 term of the U{1) gauge field, M~ 1 4>F fJll/ F ,Jjl ', where M* is an associated 
mass scale and thought to be of order the Planck scale. A change in the background 
value of 0, can be interpreted as a change of the effective coupling constant. Bekenstein 
noticed that F 2 has a non-vanishing matrix element over protons and neutrons, of order 
(10~ 3 — 10~ 2 )mjv. This matrix element acts as a source in the <fi equation of motion and 
naturally leads to the cosmological evolution of the <fi field driven by the baryon energy 
density. Thus, the change in translates into a change in a on a characteristic time scale 
comparable to the lifetime of the Universe or larger. However, the presence of a massless 
scalar field in the theory leads to the existence of an additional attractive force which 
does not respect Einstein's weak universality principle. The extremely accurate checks of 
the latter || lead to a firm lower limit on M*, M*/Mpi > 10 3 that confines possible changes 
of a to the range Aa < 1(T 10 - 1(T 9 for < z < 5 [g §. 

This range is five orders of magnitude tighter than the change Aa/a ~ 10~ 5 indicated 
in the observations of quasar absorption spectra at z = 0.5 — 3.5 and recently reported 
by Webb et al. ||. Given the potential fundamental importance of such a result, one 
should remain cautious until this result is independently verified. Nevertheless, leaving 
aside the issue regarding the reliability of the conclusions reached by Webb et al. |J, 
it is interesting to explore the possibility of constructing a dynamical model, including 
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modifications of Bekenstein's model, which could produce a large change in a in the redshift 
range z = 0.5 — 3.5 and still be consistent with the constraints on Aa/a from the results 
of high-precision limits on the violation of equivalence principle by a fifth force. It is also 
interesting to study whether the range Aa/a ~ 10~ 5 could be made consistent with the 
limits on Aa/a (J7j-[|10|, extracted from the analysis of element abundances from the Oklo 
phenomenon, a natural nuclear fission reactor that occurred about 1.8 billion years ago. We 
note that while big bang nucleosynthesis provides limits on much longer timescales, these 
limits are typically quite weak, Aa/a ~ 1CT 2 fTT[ . 



The gap of five orders of magnitude between the desirable range of 10 -5 and the bounds 
of order 10~ 10 appear to be insurmountable for any sensible modification of Bekenstein's 
theory[|. In this paper, we propose a modification of Bekenstein's idea consistent with 
experimental constraints, but relies on a large coupling between the non-baryonic dark 
matter energy density and the <fi field. 

At first, such a coupling may appear strange. Indeed, why should dark matter interact 
with the field when it is known that dark matter particles are not charged [13j and their 
electromagnetic form-factors are also tightly constrained [ftljl ? It turns out that in certain 
classes of models for dark matter, and in supersymmetric models in particular, it is natural 
to expect that would couple more strongly to dark matter particles than to baryons. It is 
easy to demonstrate this idea by a simple supersymmetrization of Bekenstein's interaction. 
In addition to the coupling of to the kinetic term, F 2 , of the gauge boson, will acquire 
an additional coupling to the kinetic term of the gaugino, M^cpx^X- If this gaugino 
constitutes a significant fraction of the stable LSP neutralino, as is often the case, the 
source of due to the energy density of dark matter turns out to be dramatically enhanced 
compared to the baryonic source, 

Dark matter source ^ 2 _ ^ 1q3 _ 1q4 _ (Ll) 

baryonic source ^baryon 

Such an enhancement factor compensates, although not entirely, for the tremendous sup- 
pression of Aa once the Eotvos-Dicke-Braginsky (EDB) limits on M* are imposed. It is 
then reasonable to study this class of models in further detail as they are numerically much 
more promising than the original Bekenstein framework. 

We note that there is another possible "strategy" to avoid the EDB constraint. One can 
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assume the existence of some extremal value ex t, in the vicinity of which only (0 — cx t) 
couples to F 2 . This type of coupling was advocated in Ref. |]15| . If the cosmological 
evolution drives close to ex t now |L5" |, i.e. at z — 0, the EDB constraints will be relaxed. 



We organize this paper as follows. In the next section we generalize the original Beken- 
stein model. In section 3, we solve the field equation for the scalar field and obtain the 



A recent publication claiming that the 10 5 change in a is realistic in this framework |12 does not 
impose the limits from Eotvos-Dicke-Braginsky experiments. 
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predictions for the change of a. In the same section, we impose experimental constraints 
and compare the results for Act with the range suggested by Webb et al. ||. In section 
4, we consider predictions for Aa in some specific models and demonstrate one model that 
passes all constraints. In section 5, we analyze the class of models with quadratic couplings 
to F 2 . Our conclusions are presented in section 6. 



2 Generalization of Bekenstein's model 

We start our analysis by formulating a generic action that includes spin-2 gravity, kinetic 
and potential terms of a modulus 0, kinetic terms for the electromagnetic field and baryons 
as well as the dark matter action, 



S = J d A X\f^g 



+ J2 ™iBm{4>))Ni + \x$>X-\M x B x (<P) X T X 

i=p,n 



+ ... (2.1) 



Throughout this paper we assume a H signature for the metric tensor. In ( |2.1| ), 

Mpi = (SkGn)" 1 / 2 = 2.4 x 10 18 GeV is the Planck mass and M* is its analogue in the 
(f) sector. Defined this way, <ft is dimensionless. Ni stands for neutrons and protons, and 
p = 7 M (<9 A1 — ieoAfj) for protons and p = j^d^ for neutrons. Here e$ is the bare charge 
which remains constant throughout the cosmological evolution (modulo the standard RG 
evolution of eo which can be neglected in our analysis). For definiteness, we assume that 
the dark matter is predominantly the non-relativistic Majorana fermion \. While it is clear 
that one can associate % with a neutralino, our approach can be easily generalized to other 
forms of cold dark matter. Ellipses stand for the omitted electron and neutrino terms, 
as well as for a number of possible interaction terms (i.e. baryon anomalous magnetic 
moments, nucleon-nucleon interactions etc.). All mass and kinetic terms are supplied with 
^-dependent factors denoted Bi(4>). In this sense, the cosmological constant term acts as a 
potential for 0. 

We shall further assume that the change of (f) over cosmological scales is small, |A0| = 
\4>(t = to) — | *C 1, where t is the present age of the universe. As such, we can expand 
all couplings around the current value of <f), which we choose to be zero, <p(t = t ) = 0, 

S a (0) = 1 + Ca0+^a0 2 

S f (0) = 1 + Cf0 + ^f0 2 (2.2) 
B m {<f>) = l + &0+i&0 2 
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S x (0) = l + C x 0+^x0 2 - 

The effective fine structure constant depends on the value of 0. As such, <p(t) and Act /a 
are directly related, 

a(4>) 



AtiB f (4)) 

— = CF0+i(^-2C|)0 2 , (2.3) 
a Z 

and we have defined Aa/a as (ckq — a(t))/a . 

The cosmological evolution of <fi follows from the scalar field equation 

Mlu^ = -M^A B' A - B'J-^F^) - (B' n m n nn + B' p m p pp) - l -B' x M x ( X T x) ■ (2-4) 

In this formula, primes denote d/dcj), and the average (...) denotes a statistical average 
over a current state of the Universe. The term with F pv F^ v can be neglected to a good 
approximation as its average is zero for photons, and its contribution mediated by the 
baryon density, J2n, P n i{i\Fp.uF ^{i), is already included in the terms proportional to B' np . 
We further note that for a Dirac fermion ip, the mass term m^ipip (and the analogous 
combination for a Majorana fermion) coincides with the trace of the -^-contribution to 
stress-energy tensor, or — 3p^. Thus, the only term, that drives in the radiation 
domination epoch when p = 3p is AqB' a (see e.g. |T7|). One can easily check that the 



change of (f> induced by this term during radiation domination will be small compared to the 
A0 developed in the subsequent matter domination epoch. Restricting eq. (|2.4|) to matter 
domination, and assuming a linearized regime ( |2.2|) , we derive the following equation of 
motion in a Robertson- Walker spacetime with scale factor a(t): 

M*(4> + 3H<P) = -p m (C m + U<P) ~ M^A (Ca + £a0) , (2-5) 
where H = a/a and ( m is defined as 

PmCm = P X Cx + Pb(XpC P + ^nCn)- (2.6) 

Here, Y p and Y n are the abundances of neutrons and protons in the Universe, including those 
bound in nuclei. We also assume that p m = p x + pb- In a more sophisticated treatment, one 
may include the contributions of electrons, the Coulomb energy stored in nuclei and other 
minor effects. As discussed in Refs. [0, |J, to good accuracy, Cm remains constant during 
the matter dominated epoch. 

If the ^-dependent energy density becomes comparable to p m or p A = M|jA, eq. (|2.5|) 
must be solved along with Einstein's equations and energy conservation as a coupled set of 
equations. However, the small cf) solutions that we are interested in imply that p<f> is small 
and (|2~5|) can be treated separately, with a{t) used as an input function. 
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3 Cosmological evolution of the fine structure con- 
stant and the EDB constraint 



The cosmological evolution of 
becomes 

<j> + 3H<p -- 



can be determined by the Q terms in eq. (|2.5|) which 



Ml 



[CmPm + CaPa] 



M? 



C Q, 



+ Ca^/ 



(3-1) 



Here p c = ?>HqM^> x is the critical density of the Universe at t — to and f2« = Pi/p c - The 
solution to this equation can be easily found [|, [TI| 0. Throughout this paper we shall 
assume that the Universe is flat and is presently dominated by non-relativistic matter and 
a cosmological constant, fl m + = 1- in this case, the time dependence of the scale factor 
is given by 



a{tf 



and eq. (|3.1|) can be integrated analytically. The first integral is given by 



6 m m a 3 



CrJ + %(sinh(26t) - 2bt) - t c 
46 



(3.2) 



(3.3) 



where b = |fi A Hq. In principle, the constant of integration t c could be kept arbitrary. 
There is, however, only one natural way of fixing it by imposing initial conditions for deep 
inside the radiation domination epoch, i.e. at t close to 0. As discussed in the previous 
section, during radiation domination, the r.h.s of ( |3.1| ) is effectively zero. This leads to a 
~ a~ 3 scaling behavior and means that any initial value of will be efficiently damped 
by the Hubble expansion over a few Hubble times. Thus, for the solution in the matter 
dominated epoch we can safely take <fi(t = 0) = or equivalently t c = 0. 

Integrating (|3.3|) gives as a function of time, 



0(t) 



AM? 



pi 



3 M3 



!% - Cm) (Wo coth(feo) - bt coth(fet)) - Cm In Sm ^!\ 
2 smn(or J 



(3.4) 



Figure 1 shows three different types of solutions for Aa /a as a function of the red-shift 
z, where 1 + z = a /a. In this plot, we have chosen ( F = 1CT 5 , £7a = 0.7 and fl m = 0.3. 
Comparing the three curves, one can see that the variation of a at high red-shifts is mostly 
determined by ( m . If (f is negative, one would need to choose negative ( m in order to get 
smaller values of a in the past. Opposite signs of (f and ( m lead to the larger values of a 
in the past. 

Given the large parameter space, (M*, (p, Cm, Ca), one could expect that it is easy 
to get Aa(z = 0.5 — 3.5) /a ~ 10~ 5 as suggested by the analysis of the quasar absorption 
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Figure 1: Three qualitatively different types of solutions for Aa(z)/a that give smaller 
values of a in the past for positive £>. They correspond to the choice of (f = 1CT 5 and (a) 
Cm = 1, Ca — (b) Cm — 1, Ca = —2 and (c) Cm = 0, Ca = 1- The interval of z, considered 
by Webb et al., 0.5 < z < 3.5 is shown by two vertical dashed lines. 



spectra by Webb et al. |J. On the other hand, it is clear that the EDB constraints 
should severely restrict the parameter space of our model. The differential acceleration of 
two elements with different A 12 and Z 12 towards a common attractor can be expressed in 
terms of Q and to = M*/2M^ (See, e.g. 0, |), 

Ag_ = ? g(A h Z x ) 
9 



g(A 2 ,z 2 ) l 



g(Ax,Zx)+g(A 2 ,Z 2 ) 



00 



7 x 10~\ F 



Z 2 /A 1 / 3 A 



A 



-z z 

A Cn+ A Qp 



■z 1 



'3_ 

A /3 



zl 

At 



(3.5) 



where Z and A represent average Z and A of the common attractor, Z = J2 ntMiZi/ J2 n%Mi. 
The terms proportional to (f correspond to the electromagnetic contribution to the total en- 
ergy of nuclei. The best constraints on long-range forces are extracted from Ag/g measured 
in experiments that compare the acceleration of light and heavy elements. The differential 
acceleration of platinum and aluminium is < 2 x 10~ 12 at the 2a level (last reference in [|| 
as quoted in ||), and the differential acceleration of the Moon (silica-dominated) and the 
Earth (iron-dominated) towards the Sun is < 0.92 x 10~ 12 ||19|| . Choosing the appropriate 
values of Z and A and retaining only the hydrogen contribution to the mass of the Sun, we 
get 



1 

CO 

1 

CO 



Cp(Cn 
Cp(Cn 



2.9 x 10- 2 CX 



< 2.5 x 10" 11 



c P + i.8x io- 2 Cf; 



< 2.5 x 10 



-ii 



Al/Pt system 
Si/Fe system 



(3.6) 
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Log 10 (C m /V(0) 



Figure 2: The (( m /y/uj, (f/V^) parameter space. The dark-shaded (green) region is 
consistent with both the EDB constraints and with a possible relative change of a at the 
10~ 5 level, as suggested by Webb et al. 0. The light shaded (blue) region is excluded by 
EDB constraints. £a is set to zero in this plot. 



These limits were also considered in a recent paper |2(J . ( n — ( p and (p enter in eqs. §T 
in different linear combinations. Thus, it is possible to extract separate limits on uj~ 1 ( p (f 
and i^ _1 Cp(Cn — Cp)- Models that have non-zero (p also have non- vanishing £ p n unless some 
intricate conspiracy of quark, gluon and photon contributions occur. Barring such possible 
cancellations, one obtains \( n , P \ ^ \(n — ( P \ ^ 10~ 3 |Cf|. Using these relations, we can 



combine the preferred range of Ref. with the constraints, imposed by eqs. (|3.6|). 

The region excluded by the EDB constraints in the (Cm/v^ Cf/V^) parameter space 
is shown by the light shaded (blue) region in Figure |2|. Here we have set Ca = . The long 
negative-sloped band that connects the upper-left and lower-right hand corners is the range 
that reproduces Aa/a = 10~ 5 in the interval 0.5 < z < 3.5. In the original Bekenstein 
model, Cm = (10 -4 to 1CT 3 )Cf and corresponds to the positive sloped band close to the 
upper-left corner ^ As one can see, the diamond-shaped intersection is deep inside the 
range excluded by the EDB experiments. Of course, this is in agreement with conclusions 
of p. Finally, the dark-shaded (green) area represents the choice of parameters that can 
reproduce Aa/a = 10~ 5 || and still be in agreement with the EDB constraints. For this 
region, Cm/v^ ~ 3 x 10~ 3 and Cf/"\A^ < 10 -3 , which points towards models in which 
couples to dark matter and the couplings to baryons and (f are suppressed. 

In addition, we must check whether or not these choices of parameters which satisfy 

2 C™ = lO -3 ^ would require rather "generous" assumptions concerning nucleon matrix elements and/or 
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Figure 3: A plot of log 10 (|Aa(z = 0.07)/Aa(> = 0.5) | ) (a) and log 10 (|Aa(« = 
0.07)/Aa(2 = 3.5)|) (b) as a function of Ca/C™ for the choice of Q\ = 0.7 and fl m = 0.3. 
The portion of the curves below the horizontal dashed line are consistent with the Oklo 
limits || and the Webb et al. suggested change in a. 

the EDB constraints are also in agreement with limits on Aa/a, derived from isotope 
abundances in the Oklo natural reactor. Typically, these limits are strong, |Aa/a;| < 
1.2 xlO -7 || and go back to z ~ 0.140. This seems to be dramatically smaller than the range 
suggested by 0. Moreover, there is no way of suppressing Aa(z < 0.14)/Aa(0.5 < z < 3.5) 
below the 10~ 2 level using our freedom in or uo, as these parameters cancel in the ratio. 

There is, however, an extra free parameter which may be used in an attempt to reconcile 
a change of 10 -5 at 0.5 < z < 3.5 and the Oklo limit. The behavior of curve (b) in 
Figure 1 suggests that (a can be used to make Aa almost flat at z < 0.2. In order to 
determine the requirements on we quantify the comparison between "Oklo change" and 
"quasar change" as follows. In the case of the Oklo constraints, in principle, one needs 
to average a(t) over the interval < to — t < 2 x 10 9 yr. Since the exact timing of Oklo 
event is known only approximately, we choose to quantify it by simply taking a at the 
half of the Oklo redshift, Aa(z = 0.07)/a. This value must be approximately two orders 
of magnitude smaller than Aa/a, suggested by Webb et al. Thus, we consider the ratio, 
\Aa{z = 0.07)/|Aa(* = 0.5) | and \ Aa{z = 0.07)/|Aa(^ = 3.5) | as a function of ( A /( m - The 
logarithms of these ratios are plotted in Figure 3. As one can see, these ratios are two funnel- 

3 The redshift, z = 0.014, corresponds to the choice fl\ = 0.7, f2 m = 0.3, ho = 0.65, and we have assumed 
that the Oklo event took place 1.8 Gyr ago. 
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like curves and it is possible to choose (a/ Cm in such a way that Aaokio/Aa quasa r < 1CT 2 . 
For z = 3.5 this can be done rather easily in the range —2.2 < (\/(m < —1.2. For z = 0.5 
one has to choose this ratio rather carefully, (a/ Cm — — 1-7, and requires a 5% -10% fine- 
tuning. We also note that this specific value of Ca/C™ is ver Y sensitive to the choice of 
and Q m and varies significantly when Q\ and Q m are varied within their current error bars. 
The use of more restrictive bounds from Oklo by |J would only worsen the fine tuning. 

The above exercise allows us to conclude that in principle a generalized Bekenstein-like 
model can yield Aa/a ~ 10~ 5 at 0.5 < z < 3.5 and still be in agreement with the EDB 
and Oklo constraints. The limits from non-universality of a fifth force could be evaded in 
models with large couplings to dark matter and small couplings to baryons and F^F^ . 
Also, the Oklo bounds could be avoided or softened if the dark matter provides a negative 
push to at later epoch. Although such a suppression of Aaokio niay happen, it would 
appear to be highly accidental. 



4 Model realizations 

It is important to note that neither the original Bekenstein model || nor its modifications 
discussed here are fully defined at the quantum level. Indeed, the Bpi^F^F^ term con- 
tains not only the bare QED Lagrangian term but also higher dimensional operators such 
as ^PFpvFpv. It is clear then that at the loop level this will create all other possible inter- 
actions such as <f) n m e ee, <f) n m q qq, etc., generally, with divergent coefficients which cannot 
be fixed from first principles. While these terms are not expected to drastically change 
the model if one makes some plausible assumptions about the cutoff in the theory, there 
is, however, the set of operators contained in B\(<j)) which are very sensitive to the cutoff 
and are very important as they can give rise to the mass of the field, an effective cos- 
mological constant, etc. Unfortunately, the present status of the underlying theory does 
not allow for a meaningful calculation of B\((j))- This problem is, of course, tightly related 
to the cosmological constant problem | 2lfl , and/or to the smallness of the mass term for 



the quintessence field. As we have nothing to add to these issues, we must assume that 
(a and £a are basically incalculable input parameters and fix A to its value implied by the 
observation of high z supernovae, the anisotropy of the cosmic microwave background and 
large scale structure formation. In what follows, we compile a list of models which predict 
certain values for the Q couplings and/or uj and confront them with the phenomenological 
constraints, discussed in the previous section. 

1 The Bekenstein model 
In this model, one initially introduces the coupling of to F 2 , Bp ((f)) = exp(— 20) or 
(f = —2. B\ can be set to a constant so that £a = 0. In the original model, u = 1, 
however, we will keep it arbitrary for now. The change of is driven by the electromagnetic 
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fraction of the baryon energy density. The coupling of <p to nucleons is given by the same 
matrix elements, ( N = m^ 1 (N\^fF^F^\N) ~ -m^ 1 {N\&(E 2 - B 2 )\N), that determine 
the contribution of a "photon cloud" to the nucleon mass. Both the Naive quark model 
and dispersion approaches give consistent estimates of these matrix elements f22fl . Using 
the results of |22|| , presumably valid to 50% accuracy, we find that ( p ~ — 0.0007Cp and 



( n ~ 0.00015Cf- Incidentally, these values almost coincide with simple extrapolations of 
the nuclear mass formula to Z = 1, 0: ( p ~ —0.0007, ( n = 0. Since (b is determined mostly 
by Cp; Cm = Cb(^b/^m) ~ — 10~ 4 Cf- As we have discussed earlier, the constraints from EDB 
experiments, as exemplified in Figure |2], do not allow a to change by more than 1 part per 
billion at red-shifts z < 3.5. 



Restricting our attention to small variations in a, we see from eq. ( p.3[ ) that 

Aft 



a 



Cf0 (4.1) 



Then evaluating eq. (|3.4|) at z = 3.5, we find that 



— = — CfU * -lO-^Cl (4.2) 
ft uo 

Note that in this model (p and ( m are of opposite sign and the final result does not depend 
on the sign of (p. Thus from eq. ( f4.2| ), we see that this model leads to larger values of a in 
the past, which is opposite the trend reported by Webb et al. Moreover, from Figure ^, we 
see that the EDB constraint requires that ( F /u < 10~ 6 and thus we see again that | Aa/a \ 
is limited to O(10~ 10 ), in agreement with the results of ||. These results however, differ 
from those of Ref. [12| in both the allowed magnitude and sign of Aa(z < 3.5) /ft. 

2 A String- dilaton-type model 



The starting point for this class of models is the action 

J d 4 xV=g~ exp(- V2(j>) (r + {d^f + A + £ mattcr ) . (4.3) 

The functions -B«(0) are easily obtained by making a conformal transformation to the Ein- 
stein frame. We find that (f = —a/2, Ca = \/2, and ( m = y/2/2. Furthermore, since there 
is only one scale in the theory, the Planck scale, we have u = 1/2. Therefore, we are able 
to obtain a definite value, Aa/a ~ —3, (over the redshift range z = 3.5 to 0. Clearly this 
is not realistic and is related to the well known problem of a massless runaway dilaton in 
string theory. Moreover, such a model is ruled out by the EDB constraints, as it predicts 
£ np ~ 1 and ( n — ( p ~ 10~ 3 . Until more can be said about the function B\((p), there is no 
useful way to use string theory to predict changes in a. 
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3 A Brans-Dicke model 



In a Brans-Dicke model is initially coupled only to the gravitational sector of the 
theory 



Sbd = — n \ d x^f^g 
IottGat 



+ ^matter (4.4) 



As in the previous example, it is easy to show that after the conformal rescaling of metric 
to the standard Einstein frame and a field redefinition for 0, the "new" field acquires 
a universal coupling ( to the mass sector of the matter fields which is given in terms of 



'■ Cm = — 1/ \J^>J + 6. One also obtains (a = — 2y 2/ {2lo + 3). As in the Bekenstein model, 
all physical results depend only on (/ y/uj. 

Due to conformal invariance of the action for the gauge field, a tree level coupling of to 
F 2 is absent yielding = 0. However, the conformal symmetry is anomalous, and at the 
one loop level a <j)F 2 term can be generated. In some sense, the couplings of to the quarks 
and leptons will be similar (apart from their magnitude) to those of the Higgs boson. It is 
then clear that a non-zero value of will be generated through the loops of charged particles 
in the same way that the Higgs-7-7 coupling is generated. For example, the coupling of 
the Higgs boson, h, to F 2 due to the top quark loop can be obtained by differentiating 
the top quark contribution to the QED /^-function, F 2 \n(A uv / (m t (l + h/v))) — > —F 2 h/v, 
where v is the Higgs vacuum expectation value. This assumes that the ultra-violet cutoff 
Auv is h- independent. In principle, there could be different possibilities with regard to the 
0-dependence of IVjjv- If one assumes that the regulator mass depends on in exactly the 
same way as an ordinary mass, then drops out of the loop amplitude, (f is not generated 
and the Brans-Dicke scalar respects the weak equivalence principle even at the one-loop 
level. A different result would arise if we postulate a 0-independent regularization Ajjv I n 
this case, a non-zero value for (p is generated and one would typically have 

Cf = ~—Cm (-- + 1 + 2 x ~) ~ -1.5 x KT 3 C„. (4.5) 



7T V 4 9 

The three terms in parenthesis correspond to the contributions from IV-bosons, charged 
leptons and quarks from second and third generations, (u and d quarks require a separate 
and quite complicated treatment. However, their main contribution is given by the charged 
pion loop and turns out to be numerically small compared to the contributions of heavy 
quarks.) The couplings of to baryons will be simply ( m , and its non- universality appears 
at the ( n — ( p ~ 10~ 3 Cm level. Thus, the bounds ( |3.6|) push the constraints on Cm 00 " 1 to the 
level of 10~ 8 or so and leave no room for a O(10~ 5 ) relative change of a at 0.5 < z < 3.5. 
The maximum allowed change is not expected to be larger than 10~ n . The result of Webb 
et al. cannot be accommodated in a Brans-Dicke model. 
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4 A Supersymmetrized Bekenstein model 



While there are many different possible supersymmetric generalizations of the original 
Bekenstein model, we consider the simplest version which begins by promoting B F (<p)(FF) 
to the rank of a superpotential: 

- J d^Bp^F^F^ -> J d A xd 2 B 1 1 B F Q>)W v yV il + (h.c.).. (4.6) 

Here <p denotes a chiral superfield, which has <fi as its bosonic component and W is the 
supersymmetric field strength. In component notation this interaction can be rewritten as 

J d 4 x (B F ((f>) [~\f, u F, u + \xd„%x] ~ ^B'MF^x) - (4-7) 

Fff, denotes the F-component of <p and ellipses stand for other terms not relevant for the 
present discussion. We see that in addition to the interaction with the gauge boson, (ft 
acquires a coupling to the gauge fermion or gaugino, x- F^ may acquire a v.e.v. which 
contributes to the supersymmetry-breaking gaugino mass. There may also be additional 
soft-breaking contributions leading to a mass term of the form \Mx T X- Performing the 
rescaling x ~^ xl \f B F (4>), we arrive at the following Lagrangian in the <t> — x sector, 



In the linearized version of the theory given by eqs. fl2.2|) , we arrive at the following 
expression for ( x , 

Cx Wx , (4-9) 

where M x = M + ( F (F^). 

Clearly, ( x can be 0(1), if ( F ~ 0(1), however its sign is not uniquely defined unless 
we make some specific assumptions about ( F , £p, (F^) and M. For example, let us take 
B F as in the original Bekenstein model so that ( F = —2. Let us further assume that 
supersymmetry breaking occurs outside the sector so that F^ = 0. In this case, M x = M 
and ( x = 2. Since the dark matter dominates the energy density of non-relativistic matter, 
we have ( m ~ ( x . Indeed, it is quite reasonable to expect that in general \( m \ ~ \( F \ which 
leads to the relation between the dark matter and baryonic sources of <fi advertized in ( |1 . 1| ) . 

The final parameter which must be specified in the model is u. In order to obtain 
consistency with the combination of EDB constraints, we must have \( yF /^Ju\ < 10~ 3 or 
uj > 4 x 10 6 . If we again assume (\ = 0, we can compute the change in the fine structure 
constant (from z = 3.5 to z = 0) 

Aa 1.2 > > . , 

= —CrCm ^ S u 4.10 

a uj 
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Type of model 


Cf 




Cm 


^9 

9 


1 Aq 1 

| a |max 

ai u.o 2; * _i o.o 


Bekenstein model 


-2 


IO" 3 


10" 4 


IO" 6 


io- 10 


"String dilaton"* 


-V2 


-V2/2 


-a/2/2 


lO- 3 


1 


Brans-Dicke model** 





—\j\jAu) + 6 


+ 6 


10 — 3 


10 11 


SUSY BM 


1 


10- 3 


1 


1CT 6 


IO" 6 


M x -driven 


IO" 4 


KT 7 - IO" 5 


0.1 


io- 12 

CO 


10~ 5 - IO" 4 



Table 1: Order of magnitude model predictions for the set of relevant couplings, non- 
universality of ^-exchange, and maximum allowed |Aa/a| at 0.5 < z < 3.5. (*) The 
tree-level form of Bi(<p) is assumed. (**) A ^-dependent cutoff is assumed. 



Note again the sign of A a predicts that a was larger in the past, although this conclusion 
could be modified if contributes to supersymmetry breaking (so that ^ 0). Also, 
because of the EDB constraint, the relative change of |Aa/a| at z in the interval 0.5 — 3.5 
would typically be at the level of IO -6 , unless some additional fine-tuning is introduced. 
(For example, if a partial cancellation between M and C,f{F^) in M x occurs, one can get 
I Cm | > | OH an d thus satisfy the constraints shown in Figure 2.) 

In contrast with the non-suspersymmetric version of the Bekenstein model, the change 
in <p from the time of the radiation domination/matter domination transition to the present 
epoch can be of order 1 or even larger. In this case, obviously, the linearized approach to 
E>i(<f)) may fail for z ^> 1. Therefore, it is impossible to determine the total change of a from 
the BBN epoch, without specifying the complete functional form for both B x (4>) and Bf(4>)- 
Nevertheless, it can be shown that if the Bi(4>) are dominated by the few first terms in the 
Taylor expansion up to z ~ 10 5 , the change of a is within the BBN bounds. Large changes 
in (p may also entail a non-negligible backreaction of the ^-dependent stress-energy tensor 
on Freidman's equations. In this case, one could get interesting effects in the expansion of 
the Universe due to the B x ((j))M x xx term, which can be interpreted at the same time as 
varying mass dark matter [23| or the potential term for 0, that has an overall factor of p m . 

An interesting consequence of models where couples to dark matter is the non- 
universality of the free fall towards an attractor dominated by dark matter, e.g. the center 
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of a galactic halo. The magnitude of differential acceleration of a system of heavy/light 
elements towards a dark matter attractor is enhanced compared to the acceleration towards 
the Sun by the factor Cm/Cp : 



kg _ Cm 
9 w 



(Cn Cp 



Zi 
A, 



Z2 
A 2 



7 x 10" 4 Cf 



'3_ 



31 

At, 



3 x 10" 



! CmCt 



(4.11) 



Depending on the ratio of Cm/Cpi this effect may be as large as 10~ 7 — 1CT 9 when the 
usual EDB constraints are satisfied. Unfortunately, this is much lower than the current 
experimental sensitivity achieved, ~ 0(1CT 3 ), in tests of differential acceleration towards 
galactic center 



5 A gaugino driven modulus 

It may happen that the modulus field that changes a is coupled primarily to the soft 
breaking parameters. Then the coupling to F 2 and baryons may only appear at the loop 
level. Let us suppose for simplicity that initially couples only to gaugino masses, 



C = E ^ 7 ^-^M,(1 + Cm»A t A 



=1,2,3 



(4.12) 



where the summation is over the three Standard Model gauge groups (color and weak indices 
are suppressed). We assume that the lightest supersymmetric particle is the neutralino x 
which is predominantly the bino. Therefore, M x ~ M\ and C, m — Cx — Cmi- 

We now consider the possibility that all couplings of 4> to standard model fields are 
induced radiatively. At the one loop level, the couplings with SU{2) and SU(3) gauge 
bosons will be generated, 

= -^Cm 2 (4.13) 

Of 

Cg = --Cm 3 (4.14) 

7T 

The calculation of these couplings is trivial: they are obtained by differentiating gluino and 
wino contributions to the corresponding beta functions over Mj. In the derivation of these 
couplings we assumed that the cutoff scale is ^-independent. 

After the breaking of the 577(2) x £7(1) gauge symmetry, (w induces a contribution to 

Cf, 

( f = sin 2 9 W ( W ~ -1.5 x 1(T 3 Ca/ 2 . (4.15) 

The coupling to baryons is mediated by (p as before or by (g or by the 4>m q qq operators, 
induced at the supersymmetric threshold. Typically, (g induces too large a coupling to 
baryons, Cn, P ~ (0.02 — 0.06)Cm 3) to be consistent with EDB limits and to = 0(1). Therefore, 
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one must require Cm 3 "C Cmi , Cm 2 • Wilson coefficients in front of supersymmetric threshold- 
induced 4>m q qq operators are expected to be at the level of 10~ 5 — 1CT 3 from (m± and (m 2 - 
This creates a coupling of (ft to nucleons at the level ( n , p ~ (10~ 6 — 10~ 4 )Cm 2 ? which is mostly 
due to a large matrix element of the m s ss operator and/or the two-loop induced G'j lu G aflu 
operators. m u uu and m d dd generate the difference between the couplings to neutron and 
proton at the level of ( n — ( p ~ (10~ 7 — 10 _5 )Ca/ 2 - Due to the small values of C Pi „, ( P — (n, 
and (p, the constraints based on the violation of the equivalence principle ( |3.6| ) do not lead 
to very restrictive bounds, 

C 2 M Juj < 10~ 4 - 1(T 2 . (4.16) 

The possibility of "choosing" (m 17 (m 2 an d uj creates sufficient freedom to satisfy EDB 
constraints and at the same time have Aa/a at 0.5 < z < 3.5 compatible with the Webb 
et al. result. Recall that 

Aa 1.2. t_ 1.8 x 10-3 

= — CfC™ ^ Cm 1 Cm 2 (4.17) 

a uj uj 

Indeed, to satisfy both it is sufficient to take 

U/ ' ±0.1; and Cm3«Cmi,Cm2- (4.18) 



uj Wuj 



To conclude this section, we combine all model predictions in Table 1. 



5 Models of an oscillating fine structure constant 



Finally, we turn to the case when all of the functions -B«(0) have a common extremum point 
</> ext . As was shown by Damour and Nordtvedt and Damour and Polyakov |Tj|, the matter 
energy density may serve as a cosmological attractor for 0, so that today its value is close 
to ext . In our approach, without loosing generality, we choose ex t = 0. The requirement 
of a common extremum is equivalent to condition that all linear couplings Q = 0. 

The cosmological evolution of is now given by the & couplings. There are two distinct 
regimes to consider: 0(t) = const at early times and an oscillating or runaway regime at 
late times. These two regimes are common for cosmological evolution of any quasi-modulus 
field, e.g. axion. The transition occurs when the Hubble rate drops below the effective 
(time-dependent) mass of <f>, 



m A = — 



Ao6v + 



Pm^rr 



6HI 

UJ 



\ 3 

^a£a + ^m£m ( —J 



(5.1) 



The sign of m 2 ^ determines if it is a runaway or oscillatory evolution. Here, we are interested 
in the oscillatory regime, and thus assume that mi is positive. The amplitude of these 
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oscillations red-shifts as a~ K , where 3/4<k<3/2. k = 3/2 occurs if ( |5.1|) is dominated by 
the first term, i.e. rigid mass, and k = 3/4 if the second matter-induced term is dominant 
]15| . Thus, the effective value of a also oscillates, at twice the frequency, 2m^, and with 



an amplitude decreasing as a 2k . In this regime, is it possible to satisfy the EDB and Oklo 
bounds and have Aa/a ~ 10 -5 at 0.5 < z < 3.5? 

If £,p is the dominant source of the couplings to baryons, the expected level of the 
violation of the equivalence principle is 

g u 

Of course, it is possible that the value of <fi today is close to zero, simply because in the 
oscillatory regime, = occurs regularly. This is, however, an accidental situation, and 
one would naturally expect now to be on the order of the amplitude of oscillations. On the 
other hand, the relative change of a is given by 

— = Uf(<P{z) - 2 now ) * (5.3) 
a Z I 

Using the relation between <p(z) and now , and plugging in the constraint from (|5.2| ), we get 

— < 10- 5 (l + z) 2K f (5.4) 

It is then clear that this can be consistent with 10~ 5 at 0.5 < z < 3.5 naturally without 
a fine-tuning of parameters when uj/^f ~ 1- The Oklo bounds can be made marginally 
consistent with |J in this scenario only for large z (close to 3.5 rather than 0.5) and for 
large k, k = 3/2. This favors models where the oscillations of <fi at z < 3.5 are driven by 
the rigid ^A-proportional mass term of 0. 



6 Conclusions 

We have shown that generalized versions of the Bekenstein's model may be consistent with 
the strong limits imposed by Eotvos-Dicke-Braginsky type of experiments and at the same 
time provide a relative change of a at the 10~ 5 -level, claimed recently in Webb et al. |J. 
The necessary flexibility in our models is achieved by the coupling of the modulus field <fi 
to the dark matter energy density and to the cosmological constant. We argue that it is 
natural to expect that the cosmological evolution of <fi will be mostly driven by these sources 
rather than by the baryon energy density. This can be seen explicitly in the simplest SUSY- 
version of the Bekenstein model, where the supersymmetric partner of the U(l) gauge field 
is the dominant non-baryonic component of dark matter. 
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In practice, it turns out that among various models where couples to F 2 , baryons 
and dark matter, only a few survive the EDB constraints and provide the O(10~ 5 ) relative 
change in a over the redshift range 0.5 < z < 3.5. In particular, we find that the models 
where is coupled initially only to £7(1) and 577(2) gaugino mass terms can easily satisfy 
both criteria. 

The bounds on Aa from the Oklo phenomenon are less dependent on the details of the 
coupling of to the matter field. Generally, they are strong enough to rule out the change 
of the fine structure constant, implied by Webb et al. In the context of the generalized 
models discussed here, the negative coupling of to the cosmological constant may be used 
to slow down its evolution and make Oklo bounds consistent with ||. This possibility, 
however, looks accidental and very fine-tuned for Act/a ~ 10~ 5 at z — 0.5. Of course, our 
treatment of all models at the loop level is plagued by the usual problem of the cosmological 
constant and the near-masslessness of the moduli field 0. This prevents us from making 
any prediction for the size of the (\ coupling constant. We also find that Q = |L5] is easier 
to reconcile with EDB constraints and Webb et al., as in this case there is an additional 
suppression of the 0- mediated force. 

Acknowledgments 

This work was supported in part by DOE grant DE-FG02-94ER-40823 at the University 
of Minnesota and NSERC. M.P. would like to thank the theory groups at McGill University 
and UQAM for their warm hospitality. We benefited from conversations with R. Cyburt, 
D. Demir, B. Fields, T. ter Veldhuis and especially M. Voloshin. 



References 

[1] P.A.M. Dirac, Nature, 139 (1937) 323. 

[2] P. Sisterna and H. Vucetich, Phys. Rev. D 41 (1990) 1034. 

[3] J. D. Bekenstein, Phys. Rev. D 25 (1982) 1527. 

[4] R.V. Eotvos, V. Pekar and E. Fekete, Ann. Phys. (Leipzig) 68 (1922) 11; 
P. G. Roll, R. Krotkov and R. H. Dicke, Annals Phys. 26 (1964) 442; 
V.B. Braginsky and V.I. Panov, Zh. Eksp. Teor. Fiz. 61 (1972) 873 [Sov. Phys. JETP 
34 (1972) 463. 

[5] M. Livio and M. Stiavelli, Ap. J. Lett. 507 (1998) L13. 
[6] J.K. Webb et al, Phys. Rev. Lett. 87 (2001) 091301. 
[7] A. I. Shlyakhter, Nature 264 (1976) 340. 



17 



T. Damour and F. Dyson, Nucl. Phys. B 480 (1996) 37. 
[9] Y. Fujii et al, Nucl. Phys. B573 (2000) 377. 



[10] S. J. Landau and H. Vucetich, |astro-ph/ 00053161 N. Chamoun, S. J. Landau and 



H. Vucetich, Phys. Lett. B 504 (2001) 1. 

[11] E. W. Kolb, M. J. Perry and T. P. Walker, Phys. Rev. D 33, 869 (1986); B. A. Campbell 
and K. A. Olive, Phys. Lett. B 345, 429 (1995) [|arXiv:hep-ph/9411272|[ ; L. Bergstrom, 
S. Iguri and H. Rubinstein, Phys. Rev. D 60, 045005 (1999) ||arXiv:astro-ph/9902157 



[12] H. B. Sandvik, J. D. Barrow and J. Magueijo, |arXiv: astro-ph / 1 075 1 2| ; J. D. Barrow 



H. B. Sandvik and J. Magueijo, |arXiv:astro-ph /0109414| . 



[13] S. Davidson, B. Campbell and D. Bailey, Phys. Rev.D43 (1991) 2314. 

[14] M. Pospelov and T. ter Veldhuis, Phys. Lett. B 480 (2000) 181 ||hep-ph/0003TJlT] 



[15] T. Damour and K. Nordtvedt, Phys. Rev. Lett. 70 (1993) 2217; 
T. Damour and K. Nordtvedt, Phys. Rev. D 48 (1993) 3436; 
T. Damour and A. M. Polyakov, Nucl. Phys. B 423 (1994) 532 ||hep-th/940i069 



[16] N. Kaloper and K. A. Olive, Astropart. Phys. 1 (1993) 185. 

[17] A. A. Tseytlin, Mod. Phys. Lett. A6 (1991) 1721; Class. Quant. Grav. 9 (1992) 979; 
Int. J. Mod. Phys. Dl (1992) 223. 

[18] R. H. Dicke, The Theoretical Significance of Experimental Relativity (Gordon and 
Breach, New York, 1965). 

[19] J.O. Dickey et al, Science 265 (1994) 482; 

J.G. Williams et al, Phys. Rev. D53 (1996) 6730; 



T. Damour, |arXiv:gr-qc/0109063 



[20] G. R. Dvali and M. Zaldarriaga, |hep-ph/0108217| . For some unknown reasons, the 



(dominant) contribution from the nuclear electromagnetic energy was not taken into 
account. 

[21] S. Weinberg, Rev. Mod. Phys. 61 (1989) 1. 

[22] J. Gasser and H. Leutwyler, Phys. Rept. 87 (1982) 77. 

[23] G. W. Anderson and S. M. Carroll, |arXiv:astro-ph/971 1288 . 



[24] G. Smith et al, Phys. Rev. Lett. 70 (1993) 123. 



18 



